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Abstract:  Compressive  sensing  (CS)  is  compared  with  conventional 
beamforming  using  horizontal  beamforming  of  at-sea,  towed-array  data. 
They  are  compared  qualitatively  using  bearing  time  records  and  quanti¬ 
tatively  using  signal-to-interference  ratio.  Qualitatively,  CS  exhibits 
lower  levels  of  background  interference  than  conventional  beamforming. 
Furthermore,  bearing  time  records  show  increasing,  but  tolerable,  levels 
of  background  interference  when  the  number  of  elements  is  decreased. 
For  the  full  array,  CS  generates  signal-to-interference  ratio  of  12  dB,  but 
conventional  beamforming  only  8  dB.  The  superiority  of  CS  over  con¬ 
ventional  beamforming  is  much  more  pronounced  with  undersampling. 
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1.  Introduction 

Beamforming  has  achieved  wide  success  in  passive  sonar  applications  such  as  the 
detection  and  localization  of  sound  sources.1’2  Plane  wave  beamforming  is  commonly 
applied  to  ship-towed  horizontal  arrays,  where  it  is  used  to  detect  targets,  determine 
target  bearing,  and  enhance  signal-to-noise  ratio.  Specifically,  it  is  a  technique  for  the 
estimation  of  the  spatial  Fourier  wavenumber  spectrum  from  measurements  of  a  spa¬ 
tially  varying  acoustic  field.  Discrete  peaks  in  the  spectrum  are  associated  with  sources 
of  sound  localized  in  space.3  However,  beamforming  is  a  conventional  l2  technique 
and  requires  that  measurements  be  properly  sampled;  i.e.,  the  array  must  have  at  least 
twice  as  many  samples  as  the  highest  wavenumber  in  order  to  prevent  aliasing.  In  this 
paper,  the  compressive  sensing  technique  is  applied  to  bearing  estimation  using  a 
towed  array.  Results  demonstrate  that  compressive  sensing  provides  better  estimates  of 
target  wavenumbers  (bearing)  while  only  requiring  approximately  one  eighth  of  the 
number  of  measured  samples  as  conventional  techniques. 

In  compressive  sensing  theory,  /] -minimization  has  proven  to  be  an  efficient 
method  for  exactly  reconstructing  the  sparsest  solution  to  underdetermined  problems 
of  linear  equations.4  Let  there  be  an  unknown  signal  of  interest,  x0  €  R^.  The  signal 
is  sparse:  it  has  only  a  few  values  that  are  significantly  larger  than  the  rest.  It  is  pos¬ 
sible  to  generate  a  vector  of  measured  data,  y  €  IM,  via  the  linear  projection 
A x=y,  where  the  sensing  matrix  A  is  chosen  such  that  it  has  rank  M  and  dimension 
M<N.  Unlike  conventional  l2  techniques,  this  problem  can  be  solved  even  though 
the  measurements  in  y  are  fewer  than  the  dimension  of  x  when  xa  is  sufficiently 
sparse — namely  that  the  number  of  nonzero  entries  of  x  is  K,  and  K  <C  M.  The 
sparsity-sensing  properties  of  this  technique  have  made  it  ideal  for  such  applications 
as  bandwidth  compression,  image  recovery,  and  signal  recovery.5’11  In  this  paper  an 
alternating  direction  algorithm  is  used  to  solve  this  convex  optimization  problem.7 
Specifically,  the  YALL1  algorithm  was  selected  after  comparison  with  other  cur¬ 
rently  available  algorithms,  and  found  to  produce  repeatable  inversion  results  for 
beamforming  using  compressive  sensing.8  The  interface  of  YALL1  enabled  simpler 
implementation  with  existing  code. 
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2.  Beamforming  using  compressive  sensing 

An  algorithm  is  proposed  for  passive  bearing  estimation  where  the  number  of  sources 
and  their  source  signals  are  considered  unknown  in  the  presence  of  noise.  An  unknown 
but  limited  number  of  targets  is  assumed  to  be  distributed  in  bearing  (look  angle)  with 
unknown  bearings  and  magnitudes.  The  receiver  consists  of  a  linear  array  (either  hori¬ 
zontal  or  vertical)  that  has  arbitrary  inter-element  spacing.  Let  r  be  the  vector  of  re¬ 
ceiver  positions.  The  classical  discretized  equation  for  a  beamformer  is 

Amn  =  wm  exp [ik  sin ((j>n)rm\  (1) 

where  A  is  the  phase  and  amplitude  weights  calculated  for  any  frequency  f=kc!2n, 
and  c  is  the  soundspeed  of  the  medium.  The  look  directions,  </>„,  range  from  0  to  n 
radians  with  nil  being  dehned  as  broadside  to  the  array.  Although  not  necessary,  the 
array  may  be  shaded  by  w  (e.g.,  harming,  hamming,  or  box  car).  For  an  M  element 
array,  N  look  directions  are  assumed  in  the  beam  domain.  Here  the  seemingly  under- 
determined  problem  where  N  >  M  is  considered.  The  sparsity  solution  x  comprises  the 
angle  and  amplitude  information  of  the  A  <  M  targets  present  and  can  be  estimated  by 
the  following  l\  minimization  problem: 

x  =  min  ||  x  || y  subject  to  Ax  =  y ,  (2) 

where  y  is  a  vector  of  measurements  y  =  \y\,  y%,—,  >’„,]7  and  x  is  the  solution  of  the 
unknown  x. 

2. 1  Broadband  estimation  of  an  unknown  source  signed 

In  passive  acoustic  problems,  the  source  function  of  the  target  of  interest  is  often  con¬ 
sidered  unknown.  Therefore,  a  coherent  broadband  A  matrix  cannot  be  constructed  to 
search  for  the  sparse  look  direction  vector  x.  The  approached  used  here  is  to  solve 
each  frequency  independently,  and  to  average  incoherently  the  solutions  producing  a 
final  beamformer  output  of  x  =  1  |*(  fj)  |  f°r  all  J  frequencies. 

2.2  Stability  of  A 

The  restricted  isometry  property  (RIP)  relates  to  the  orthonormality  of  the  columns  of 
the  A  matrix  for  sparse  recovery.9  The  RIP  condition  is  a  necessary,  though  not  suffi¬ 
cient,  condition  for  signal  recovery.  The  ^-restricted  isometry  constant  is  defined  as  the 
smallest  8k  value  that  satisfies 

(1  -  4)  \\b\\l  <  \\Ab\\l  <  (1  +  4)  ||*|||  (3) 

for  every  possible  /r-sparsc  vector  b.  Additionally,  the  k,  Id -restricted  orthogonality 
constant  ©k.k'  is  defined  as  the  smallest  value  that  satisfies 

\(Ab,Ab')\<©KK,\\b\\2\\b'\\2  (4) 

for  all  b  and  //  such  that  they  are  A-sparse  and  Ab -sparse,  respectively,  and  have  disjoint 
supports.  For  the  data  analysis  presented  in  this  paper,  the  RIP  condition  always  satis¬ 
fied  the  inversion-without-noise  criteria,  namely,  c>  i  .75^  <  \/2  -  I ,  where  k=  4  sparse 
targets  were  present.10  The  recovery  of  sparse  signals  in  the  presence  of  bounded  error 
is  possible  if  <5 1 .25*  +  ©k.x.isk  <  1-  This  condition  was  almost  met;  a  value  of  1.05  is 
obtained  for  the  fully  populated  array  case.  However,  numerical  inversions  remained 
stable  to  noise  without  use  of  the  l2  constraint  (i.e.,  noisy  compressive  sensing11). 

Most  linear  arrays  have  inter-element  spacing  designed  for  a  specific  fre¬ 
quency.  When  such  an  array  is  used  at  a  lower-than-design  frequency,  A  becomes 
ill-conditioned.  Compressive  sensing  inversions  are  seeded  with  an  l2  estimate,  therefore 
even  in  the  case  of  an  arbitrarily  spaced  array,  the  rank  of  A  may  be  insufficient, 
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leading  to  12  inversion  instabilities  and  a  poorly  defined  starting  point.  In  order  to 
improve  the  rank  of  A,  while  retaining  potentially  valuable  data  in  y,  the  use  of  a  sin¬ 
gular  value  decomposition  (SVD)  to  intelligently  precondition  the  basis  set  and  data  is 
used  here.  The  SVD  of  A  (where  11  indicates  the  conjugate  transpose)  produces: 

USV^x  =  y.  (5) 

The  diagonal  of  S  contains  the  singular  values  of  A.  For  the  data  analysis  presented  in 
this  paper,  all  singular  values  <  1/100  of  the  maximum  value  are  set  to  zero.  Thus,  the 
size  of  the  A  matrix  will  be  reduced  from  N  x  M  to  N  x  P,  where  P  is  the  number  of 
retained  singular  values  of  S.  After  reducing  the  U  and  S  matrices  to  their  non-zero 
columns  and  rows  (denoted  by  U  and  S,  respectively): 

\}"\]SVn x  =  l]Uy,  or  equivalently  Ax  =  y,  (6) 

x  =  arg  ||  x  || ! ,  subject  to  Ax  =  y.  (7) 

Thus,  the  A  matrix  is  ensured  to  have  a  stable  condition  number  for  inversion  and  all 
of  the  (redundant)  data  are  maximally  utilized.  Additionally,  reducing  the  size  of  A 
and  y  also  has  the  effect  of  decreasing  the  computational  time  required  for  inversion 
results.  The  SVD  can  also  generate  orthonormality,  and  thereby  improve  the  RIP  con¬ 
dition,  though,  in  this  case,  the  beam  matrix  A  is  normal  as  defined. 

3.  Underwater  acoustic  results  from  the  BASE07  experiment 

Beamforming  using  compressive  sensing  is  demonstrated  on  at-sea  data  measured  dur¬ 
ing  the  BASE07  experiment.12’13  During  this  experiment,  the  Five  Octave  Research 
Array  (FORA)14  was  towed  behind  a  ship  near  the  Malta  Plateau.  The  FORA  consists 
of  several  nested  arrays  of  which,  for  the  purpose  of  this  paper,  the  mid  frequency  ele¬ 
ment  spacing  array  is  considered  (design  frequency  1000  FIz,  0.75  m  spacing,  M=64, 
f =  8000  Hz).  During  Julian  Day  135  at  approximately  13:15  UTC,  the  towed  array 
received  a  loud,  broadband  FM  broadcast  by  another  ship.  A  60  s  spectrogram  of  the 
first  element  of  the  FORA  array  is  shown  in  the  upper  left  panel  of  Fig.  1,  which 
depicts  the  FM  arrival  at  approximately  14  s.  Though  not  apparent  in  the  background 
noise  of  the  spectrogram,  at  least  four  other  ships  may  be  present. 

3.1  Conventional  beamforming  results 

Conventional  beamforming  results  are  shown  in  the  upper  right  panel  of  Fig.  1.  Beam- 
former  output  is  shown  as  a  function  of  time  and  look  angle.  The  60  s  time  series  is  di¬ 
vided  into  time  steps  of  2048  points  (0.256  s)  with  no  overlap.  At  each  time  step, 
N  =  M  =  64  number  of  beams  were  calculated  and  applied  to  the  measured  data.  The 
final  beamformer  output  is  the  incoherent  average  from  800—1000  Hz  (df=  3.9  Hz). 
The  loud  FM  sweep  can  be  seen  at  approximately  14.15  s  and  at  -31°.  Additionally, 
there  appears  to  be  at  least  4  sources  of  opportunity  made  at  bearing  angles  -61°, 
—44°,  7°,  and  32°.  One  of  the  three  sources  (7°)  clearly  changes  bearing  by  the  end  of 
the  minute.  The  signal  strength  of  the  FM  sweep  is  approximately  30  dB  stronger  than 
the  sources,  thus  the  dynamic  range  of  the  plot  is  set  to  60  dB. 

3.2  Compressive  sensing  beamforming  results 

Beamforming  results  using  compressive  sensing  are  shown  in  the  lower  right  panel  of 
Fig.  1.  The  same  time  steps  are  used;  however,  at  each  time  step,  N=  128  >  M  number 
of  beams  were  calculated  and  applied  to  the  measured  data.  Conventionally,  this  is  a 
highly  underdetermined  problem;  however  the  compressive  sensing  assumption  of 
sparseness  and  use  of  l\  minimization  requires  significantly  less  information  to  achieve 
a  unique  and  meaningful  solution  (2 N  vs  K  log  N).  Compressive  sensing  can  yield 
greater  angular  resolution  than  is  conventionally  possible  for  a  given  array  of  M 
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Fig.  1.  (Color  online)  BASE07  processed  results.  The  spectrogram  of  the  first  element  of  the  towed  FORA 
array  shows  a  loud  FM  arrival  at  approximately  15  s.  The  right  most  panels  show  beamformer  outputs  incoher¬ 
ently  averaged  from  800  to  1000  FIz.  The  upper  panel  depicts  conventional  beamforming  results  and  the  bottom 
shows  compressive  sensing  results  both  over  60  dB  dynamic  range.  Compressive  sensing  yields  smaller  angular 
resolution,  reduced  angular  spillage  of  the  strong  FM  sweep,  and  tighter  endfire  energy.  A  time  slice  of  the 
beamformer  outputs  at  5  s  is  depicted  in  the  bottom  left  panel.  Comparing  the  peaks  and  valleys  of  the  compres¬ 
sive  sensing  (lower)  and  conventional  (upper)  results  show  that  compressive  sensing  has  superior  interference 
and  noise  suppression  (note  conventional  results  raised  by  +9  dB  to  align  peak). 


elements.  This  is  seen  in  the  beam-former  outputs  in  the  right  two  panels  in  Fig.  1:  the 
ship  lines  at  about  7°  and  32°  are  more  narrow  and  focused  in  the  lower  right  hand 
panel  than  the  upper  right  hand  panel.  Compared  with  the  conventional  beamformer, 
the  loud  FM  is  more  tightly  focused  in  angle  with  less  leakage  into  other  beams.  Addi¬ 
tionally,  the  tow  ship  noise  measured  at  endfire  appears  more  tightly  contained  at 
higher  angles  (the  endfire  lobe  being  effectively  smaller).  Note  that  the  compressive 
sensing  beamforming  results  yield  an  accurate  total  receive  level,  and  it  is  possible  to 
“listen”  to  a  bearing  angle  by  taking  an  inverse  Fourier  transform  across  frequencies. 

Noise  (distributed  in  angle)  and  interference  are  significantly  suppressed  with 
compressive  sensing  beamforming.  The  bottom  left  hand  panel  of  Fig.  1  depicts  the 
beamformer  output  at  5  s  (dashed  lines  on  the  right  panels).  The  conventional  beam¬ 
forming  results  were  offset  by  9  dB  to  align  the  peak  at  7.3°.  Comparing  peaks  to  val¬ 
leys,  compressive  sensing  provides  a  greater  main  to  interference  (and  noise)  ratio. 
Conceptually,  this  is  possible  because  the  total  energy  of  the  input  signal  is  not 
required  to  be  preserved  with  l\  beamforming.  The  improvement  is  quantified  for  the 
target  at  bearing  32°  using  the  definition  of  signal-to-interference  ratio  (SIR), 


SIR  =  201og10 


max[|v(30  <  <j>  <  35) |] 
std[|.\'(25  <  </>  <  30,35  <  <£  <  40)  |] 


(8) 


J.  Acoust.  Soc.  Am.  130  (4),  October  201 1  G.  F.  Edelmann  and  C.  F.  Gaumond:  Beamforming  using  compressive  sensing  EL235 


Downloaded  16  Nov  201 1  to  74.96.254.145.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


G.  F.  Edelmann  and  C.  F.  Gaumond:  JASA  Express  Letters 


[DOI:  10.1121/1.3632046] 


Published  Online  09  September  201 1 


Fig.  2.  (Color  online)  Bearing  versus  time  results  for  the  first  25  s  of  data.  The  total  number  of  array  elements 
(not  aperture)  is  reduced  toward  the  right.  Despite  the  array  being  reduced  from  the  original  64  elements  (half 
wavelength  spacing)  to  as  little  as  6  randomly  selected  elements,  compressive  beamforming  is  able  to  detect  tar¬ 
gets  of  interest. 


The  SIR  produced  with  compressive  sensing  is  approximately  12  dB,  while  the  SIR 
with  conventional  processing  is  approximately  8  dB.  Thus,  the  gain  generated  through 
compressive  sensing  is  approximately  4  dB.  Note  also  that  at  the  peak  (FM  sweep), 
compressive  sensing  produced  a  9  dB  gain.  Thus,  with  sufficient  signal-to-noise  ratio, 
compressive  sensing  can  strongly  outperform  conventional  processing  results. 

Compressive  sensing  is  also  able  to  reconstruct  a  sparse  signal  with  fewer 
measurements  than  would  be  expected  from  traditional  sampling  theory.6  Figure  2 
shows  the  performance  degradation  due  to  reduced  numbers  of  towed  array  hydro¬ 
phones.  Each  panel  depicts  the  first  25  s  of  broadband  compressive  beamforming  with 
different  numbers  of  randomly  chosen  hydrophones.  Note  that  the  first  and  last  hydro¬ 
phones  are  always  included  to  maintain  the  array  aperture.  The  first  panel  shows 
results  using  32  out  of  64  possible  hydrophones,  with  the  subsequent  panels  depicting 
16,  11,  8,  and  6  hydrophones,  respectively.  The  results  show  time-bearing  estimates 
without  aliasing  even  when  the  array  is  reduced  to  8  hydrophones.  For  the  cases 
shown  in  Fig.  1,  the  values  of  SIR  are  12,  9,  10,  11,  and  12  dB  respectively.  Note  that 
the  color  scale  ranges  from  85  dB  to  145  dB  and  shows  the  relative  amplitudes  of  the 
various  components.  This  large  dynamic  range  does  not  clearly  show  SIR  with  fewer 
hydrophones.  Not  shown  are  the  results  from  conventional  beamforming:  This  method 
predictably  fails  to  isolate  the  line  at  32°  with  SIR  values  of  1,  1,  1,  0,  and  0  dB, 
respectively.  Note  that  if  the  elements  of  the  array  were  to  be  equally  spaced,  then  the 
spacing  would  have  equivalently  produced  an  array  of  approximately  1,  21,  31,  41, 
and  52  spacing,  respectively.  Compressive  sensing  inversion  became  unstable  when 
using  the  6  element  array,  therefore  reducing  the  array  further  became  to  intractable. 
Also  of  interest,  the  32  element  results  are  slightly  more  superior  than  the  64  element 
results.  It  is  hypothesized  that  this  is  due  to  the  equal  spacing  of  the  full  64  hydro¬ 
phone  array;  compressive  sensing  works  best  with  an  array  that  has  been  sampled  uni¬ 
formly  at  random  so  that  a  robust  uncertainty  principle  prevails.6 

4.  Conclusion 

Compressive  sensing  is  a  novel  approach  useful  to  exactly  reconstruct  a  sparse  signal 
with  fewer  measurements  than  defined  by  traditional  sampling  theory.  This  paper 
presents  the  application  of  compressive  sensing  to  beamforming  of  measured  under¬ 
water  acoustic  data  that  is  sufficiently  sparse  in  bearing  angle.  The  implemented 
algorithms  have  the  ability  to  improve  visibility  of  lines  in  time-bearing  displays.  Com¬ 
pared  with  conventional  beamforming,  compressive  sensing  shows  finer  angular  resolu¬ 
tion  and  greater  interference  suppression  when  applied  to  measured  at-sea  data. 
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Furthermore,  the  technique  is  demonstrated  to  work  on  a  sub-sampled  array  with  no 
loss  in  SIR  performance,  even  with  a  number  of  elements  that  would  equivalently  pro¬ 
duce  a  three-wavelength  spaced  array  (where  conventional  methods  require  half  wave¬ 
length).  Therefore,  the  method  allows  the  use  of  arrays  with  dead  hydrophones,  or 
allows  the  possibility  of  acoustic  arrays  to  be  manufactured  with  fewer  elements. 
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